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(57) ABSTRACT 

A method and model of predicting particle-surface interac- 
tions with a surface, such as the surface of a spacecraft. The 
method includes the steps of: determining a trajectory path of 
a plurality of moving particles; predicting whether any of the 
moving particles will intersect a surface; predicting whether 
any of the particles will be captured by the surface and/or; 
predicting a reflected trajectory and velocity of particles 
reflected from the surface. 
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PARTICLE-SURFACE INTERACTION 
MODEL AND METHOD OF DETERMINING 
PARTICLE-SURFACE INTERACTIONS 

ORIGIN OF INVENTION 5 

The invention described herein was made by an employee 
of the United States Government, and may be manufactured 
and used by or for the Government for governmental pur- 
poses without the payment of any royalties thereon or there- 10 
for. 

RELATED APPLICATION 

This application is based upon prior filed provisional 15 
patent application Ser. No. 61/061 ,033 filed Jun. 12, 2008, the 
entire subject matter of which are incorporated by reference 
herein. 

BACKGROUND OF THE INVENTION 20 

The present invention is related to the modeling of the 
interactions between particles and surfaces. 

Modeling of the transport of particles around spacecraft in 
orbit is a poorly developed science. No current models 25 
account for the way in which real particles bounce off of 
surfaces when they impact. 

The process by which particle contamination moves 
throughout a system may be divided into three phases. The 
first phase is generation and release from a surface, the second 30 
phase is transport through air or vacuum, and the third phase 
is reflection or deposition on a surface. The first two phases 
have been the focus of prior art, such as the studies referenced 
in E. N. Borson, “A model for particle redistribution during 
spacecraft launch,” Proceedings of SPIE Vol. 4774, (2002), 35 
pp 67-78, the entire disclosure of which is incorporated by 
reference herein. 

With respect to the generation and release of particles from 
a surface, many particles are generated by wear between 
moving surfaces while other particles are “pre-existing” on a 40 
surface due to contamination built up over the life of the 
surface. 

Most particle generation on spacecraft occurs during 
ground processing. Particles shed from personnel, such as 
garment fibers and skin cells, have a low mass to drag ratio 45 
and are transported throughout the clean room by air flow. 
Abrasion of materials, such as drilling, grinding, and cutting, 
produces particles with a high mass to drag ratio, and these 
particles settle close to the point of generation. 

Deployable systems such as antenna booms, solar arrays, 50 
and covers create materials after spacecraft are in orbit. 
Actuators create particles at various velocities: pyrotechnical 
devices release high speed particles, spring cutters release 
medium speed particles, and sealed wax actuators release 
slow particles. 55 

Particles may be released from surfaces by vibration or 
other external forces. 

With respect to the transport of particles, drag is the viscous 
friction effect exerted on a moving particle by the medium it 
is in. The result is that large particles tend to stray from the 60 
flow direction because of momentum but small particles fol- 
low the flow direction and are less likely to impact surfaces. 
The equations of motion of particles suspended in air may be 
found in A. Busnaina and X. Zhu, “Submicron Particle 
Motion and Deposition in a CVD Chamber,” Institute of 65 
Environmental Sciences Proceedings Volume 2 (1992), pp 
1-9, the entire disclosure of which is incorporated by refer- 


2 

ence herein. As air density drops, the drag also drops, until in 
vacuum the particles are affected only by gravity. More infor- 
mation regarding particle motion in a vacuum can be found in 
A. L. Lee, “Particle Dispersion Around a Spacecraft”, Ameri- 
can Institute of Aeronautics and Astronautics, Aerospace Sci- 
ence Meeting, 21 st , Reno, Nev., Jan. 10-13, 1983, the entire 
disclosure of which is incorporated by reference herein. 

Regarding the reflection or deposition on a surface, when a 
moving particle impacts a surface, energy is lost through the 
compression of the particle and the surface. Some theoretical 
treatment of this phenomenon exists in the literature such as 
C. Thornton et al., “A theoretical model for the stick/bounce 
behavior of adhesive, elastic -plastic spheres”, Powder Tech- 
nology 99 (1998), pp 154-162 the entire disclosure of which 
is incorporated herein by reference. Kinetic energy is also 
converted between translational and rotational modes as the 
particle tumbles. The tendency to tumble is a function of 
surface roughness, particle shape, and impact angle. 

Whether a particle impacting a surface will stick is deter- 
mined by whether the potential energy stored in compression 
is great enough to overcome the forces holding the particle to 
the surface. Particles on a surface are held to the surface by 
several forces, but the large forces are electrostatic attraction, 
surface energy, and surface tension. 

If the particle bounces, the direction of the bounce may be 
affected by rotational momentum, surface roughness, and 
particle shape. 

The prior art is lacking in a model for this last phase of 
particle transport. The present invention provides a compre- 
hensive approach for describing and modeling the interac- 
tions between a particle and a surface. 

SUMMARY OF THE INVENTION 

The present invention is directed to a method of predicting 
particle- surface interactions with a surface such as a space- 
craft. The method includes the steps of: determining a trajec- 
tory path of a plurality of moving particles; predicting 
whether any of the moving particles will intersect a surface; 
predicting whether any of the particles will be captured by the 
surface and/or; predicting a reflected trajectory and velocity 
of particles reflected from the surface. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The features of the invention believed to be novel and the 
elements characteristic of the invention are set forth with 
particularity in the appended claims. The Figures are for 
illustration purposes only and are not drawn to scale. The 
invention itself, however, both as to organization and method 
of operation, may best be understood by reference to the 
detailed description which follows taken in conjunction with 
the accompanying drawings in which: 

FIG. 1 is a diagram of a particle — surface interaction. 

FIGS. 2 through 7 are flow charts illustrating the method 
according to the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

Referring to the Figures in more detail and particularly 
referring to FIG. 1, there is shown a diagram of a particle 
surface interaction. Arrow 10 denotes a vector of an incoming 
particle having incoming velocity v 0 . Vector 10 strikes sur- 
face 12 and is reflected from the surface as shown by arrow 14 
which denotes a vector of the reflected particle having 
reflected or outgoing velocity v x . Arrow 16 denotes a vector 
of the specular reflection of the particle. The specular reflec- 
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tion means that the angle of reflection of the reflected particle 
(as measured between the surface normal 20 and the reflected 
vector) (3 equals the angle of incidence (as measured between 
the surface normal 20 and the incoming vector) a. The dif- 
ference between the actual reflected vector 14 and the specu- 
lar reflection vector 16 is a diffuseness vector 18. The angle 
between the actual reflected vector 14 and the specular vector 
is 0. 

According to one preferred embodiment of the invention, 
there is a model for particle-surface interactions which 
includes a coefficient of restitution parameter, a diffuseness 
parameter, a coefficient of momentum transfer parameter, and 
a sticking velocity parameter. 

The coefficient of restitution (COR) is the fraction of trans- 
lational momentum retained by the particle after interacting 
with the surface. In a perfectly inelastic collision, all of the 
energy is retained by the particle and the COR is 1. The 
highest measured COR for materials is around 0.9 for glass 
beads on a glass plate. COR is a material property and applies 
to both the particle and the surface; the effective COR is the 
product of the individual particle and surface COR values. 

The coefficient of momentum transfer (COM) is a measure 
of the conversion between rotational and translational modes 
during a particle-surface interaction. COM is primarily a 
function of the particle moments of inertia. The fraction of 
energy converted on any given bounce is a random number 
with a Gaussian distribution about 0 (negative being taken as 
conversion to translational momentum, positive as conver- 
sion to rotational momentum). The COM is the standard 
deviation of the aforementioned Gaussian distribution. Nei- 
ther momentum is allowed to become negative. The conver- 
sion is applied after the energy reduction due to the COR. 

The diffuseness parameter, also known as diffusivity, is a 
measure of the randomness of the bounce direction. Alter- 
nately, the complement (1 -diffusivity) is the specularity of the 
bounce. For a perfectly specular bounce such that diffusiv- 
ity=0, the angle between the incoming direction and surface 
normal equals the angle between the reflecting direction and 
surface normal, and the reflecting vector is in the plane of the 
incoming vector and surface normal. A perfectly diffuse 
bounce has the reflecting vector distributed in a Lambertian 
fashion, with the probability of any direction being propor- 
tional to the cosine of the angle between the surface normal 
and that direction. The model of the present invention uses the 
diffusivity as an interpolation fraction between the specular 
direction and a random, Lambertian direction. Other distri- 
butions besides Lambertian may be used to model surfaces 
with anisotropic diffuseness. Diffusivity is a function of the 
particle shape and surface irregularities. The two diffusivities 
must be combined in a manner that meets certain criteria: 
when one value is 0 (perfectly specular), the result should 
equal the other value of diffusivity; the maximum value 
should be 1 ; and when both values are intermediate, the result 
should be larger than the higher of the values, because the 
higher of the values would be used if the smaller were zero. 
The method of combination that meets these criteria is to use 
the complement of the combined specularity, where the com- 
bined specularity is the product of the particle and surface 
specularities. 

An additional feature of the model is the assumption that 
conversion between rotational and translational modes 
affects the direction of travel. This is implemented by using 
the ratio of the converted momentum to the total momentum 
as an additional diffuseness term that is combined according 
to the multiplication of specularities. The result is that no 
change in diffuseness occurs when the momentum conver- 
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sion is zero, but transformation of all of the momentum would 
result in total randomization of the direction. 

The sticking velocity is the velocity below which the par- 
ticle does not possess enough kinetic energy to escape from 
5 the forces holding it to the surface. Both the particle and the 
surface contribute to the binding energy. Because kinetic 
energy is proportional to velocity squared, the effective stick- 
ing velocity is given by the square root of the sum of the 
squared velocities. The algorithm applies this test after energy 
to is converted between translational and rotational modes. If 
the particle does not possess the velocity required to escape, 
but does possess rotational velocity, the condition may be 
interpreted as rolling on the surface. 

In operation of the particle — surface interaction model, the 
1 5 sticking velocity parameter is determined to see if the particle 
sticks to the surface. If the particle does not stick to the surface 
then the coefficient of restitution parameter, diffuseness 
parameter and coefficient of momentum transfer parameter 
are determined to result in the direction that the particle 
20 bounces from the surface. 

As will be apparent to those skilled in the art of numerical 
simulation, the surface interaction model may be easily added 
to any particle transport simulation wherein the influence of 
surface bounces on the particle motion is non-negligible. 

25 According to one preferred embodiment of the invention, 
there is a method for determining particle — surface interac- 
tions. The algorithm here gives the example for motion of 
particles in a vacuum that are sufficiently separated that par- 
ticle-particle collisions are negligible. 

30 Referring to FIGS. 2 to 4, and particularly FIG. 2, the 
method according to the present invention is illustrated. The 
method includes first in step 30 inputting properties for the 
particles and surface followed by a loop in which a specified 
number of particles are tracked from generation through cap- 
35 ture on a surface in step 32, and then determining particle 
distribution statistics from the results of the tracked particles 
in step 38. 

The tracking loop of tracking particles through capture on 
a surface in step 32 includes initializing each particle to a 
40 position and velocity in step 34, using rules specified by the 
user, and then moving the particle and interacting with sur- 
faces until the particle velocity is set to zero in step 36 or the 
particle exits the simulation domain (spatial or temporal 
extent). 

45 Referring now to FIG. 3, moving the particle (step 36 from 
FIG. 2) includes calculating the path of the particle through 
space in step 40, checking for intersection of the particle path 
with surfaces in step 42, and if an intersection occurs with a 
surface. If there is no intersection with a surface, the calcula- 
50 tion is done as indicated in step 44 because there will be no 
particle — surface intersection. If there is an intersection of the 
particle with a surface, then the method continues. After it is 
determined that there is a particle — surface intersection in 
step 42, the next step 46 is checking to see if the particle is 
55 captured by the surface. If the particle is captured by the 
surface, the modeling is done as indicated by step 48 because 
there will be no further particle or surface activity. If the 
particle is not captured by the surface, then the next step 50 in 
the method is determining the new direction and velocity of 
60 the particle after the surface intersection. 

Referring now to FIG. 4, further steps are described with 
respect to calculating the new direction and velocity of the 
particle (step 50 in FIG. 3). Initially, the coefficient of resti- 
tution is determined in step 52 or may occur prior to step 30 
65 and be provided to step 52 as a constant value. More details 
will be provided hereafter on the method for determining the 
COR. 
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Then, the momentum loss is determined from the COR in 
step 54. Surface COR is combined with particle COR by 
taking the product of the two. The particle COR may be fixed, 
or may be calculated from material properties and the velocity 
at impact. 5 

In the following equations, “y” is understood to be in the 
vertical direction and “x” is understood to be in the horizontal 
direction. 

In the following equations 1 through 17, the following 
symbols have the noted meanings: 

v i , v xl =reflected velocity in the y and x directions, respec- 
tively 

v v0 , v x0 =incoming velocity in the y and x directions, 
respectively 15 

a=acceleration 
t=time 

y ^reflected height of particle 
y 0 =initial height of particle 

x ^reflected distance in x direction 20 

x 0 =initial distance in x direction 
v x0 =the initial velocity in the x direction 
v=total velocity 

u =unit vector 
v x =velocity in x direction 
y^velocity in y direction 
K=translational kinetic energy 
U=rotational kinetic energy 

m=mass of particle 30 

I=moment of inertia 
co ^rotational velocity after reflection 
oo 0 =rotational velocity before reflection 
COR=coefficient of restitution 

COM=coefficient of momentum transfer 35 

v =vector representing reflected particle 

v specula =\QcXor representing reflected particle in specular 
direction 

v lambertiarT^^ 01 representing reflected particle in a ran- 40 
dom direction 

©=angle between reflected particle vector and specular 
vector 

d=diffuseness 

The pre-impact kinetic energy is determined in step 56 as 45 
one half particle mass multiplied by the particle velocity 
(relative to the surface) squared: 


*o = 2 mv o 


Equation (1) 50 


This is multiplied by the square of the effective COR to 
determine the post impact kinetic energy in step 58: 


Ki = -m(COR z )v z . 


Equation (2) 


The coefficient of momentum transfer (COM) is deter- 
mined in step 60, or may occur prior to step 30 and be 
provided to step 60 as a constant value. More details will be 
provided hereafter on the methodology for determining 
COM. 65 

If the particle has a COM, the change in rotational energy 
is determined in step 62 as the product of the post impact 


kinetic energy and a random variable with a Gaussian distri- 
bution centered at zero with a standard deviation equal to the 
COM: 


A/(O 2 =[Gaussian(0, COM)] [ivj Equation (3) 

If the change in momentum would decrease the particle 
rotational momentum below zero, the change in momentum 
is set equal to the particle rotational momentum, indicating 
that no rotation will be present after the bounce. If the change 
in momentum would reduce the kinetic energy below a fixed 
percentage of the post impact value, the change is reduced to 
ensure that some translational energy remains; if not, the 
particle could spin in place, a non-physical result. 

The post impact kinetic energy and the particle rotational 
momentum are updated in step 64 with the change in momen- 
tum: 


K 2 =K l -AIoy 2 Equation (4) 

The post impact velocity is determined in step 66 from the 
square root of two divided by particle mass multiplied by the 
kinetic energy: 



Equation (5) 


The total energy is determined in step 68 as the sum of the 
particle rotational and translational energy. If this value is 
greater than zero, the diffuseness of the bounce is determined 
as the complement of one minus the fraction change in 
momentum divided by total momentum, multiplied by one 
minus particle diffuseness, multiplied by one minus the sur- 
face diffuseness: 


1 -d-Effective-W - ^com] [ 1 - d Particle ^ [ 1 - d Sur f ace \ Equation (6) 

The determination of the diffuseness is discussed in more 
detail hereafter. 

If the diffuseness is greater than zero, a random unit length 
direction vector is generated in step 70. The length of this 
vector is multiplied by the diffuseness. The specular reflec- 
tion direction in step 72, which was determined during the 
step of checking for intersection with surfaces in step 42 in 
FIG. 3, is multiplied by one minus the diffuseness. These two 
vectors are added together to get the post impact reflection 
direction in step 74 which is multiplied by the post impact 
velocity scalar to get the reflection velocity vector in step 76: 

Ti =T([t -dj u Specular +[d] ~u Random ) Equation (7) 

If the surface is moving, the reflection velocity vector is 
decomposed into a component normal to the surface and a 
component in the plane of the surface. The surface friction is 
used to update the in-plane component by adding the element 
in-plane motion multiplied by friction to the reflection in- 
plane component multiplied by one minus the friction. The 
surface normal component of motion is added to the reflection 
normal component. The velocity vector is re-composed by 
adding these two resultant vectors. 

A special condition exists when the surface motion is away 
from the particle; the reflection velocity may be towards the 
surface rather than away. This is interpreted as a soft landing, 
and the post impact velocity is set to zero. 

At this point, the model is complete in that the direction and 
velocity of the reflected particle has been determined after the 
particle — surface interaction. 

Referring now to FIG. 5, the methodology for obtaining the 
COR will be described. The methodology begins in step 80 
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with bouncing a plurality of particles on a plate (one at a time) 
under the influence of gravity and the location and times of 
impact are observed, for example with a pair or high speed 
cameras. 

From Newton’ s laws of motion, assuming no air drag (per- 5 
forming the test in a vacuum chamber), there is a ballistic 
trajectory: 


v n = v Y o + at Equation (8) 10 

yi =yo + v YO t + ^at 2 
xi =x 0 +v X ot 

15 

In the next step 82, the height of the bounce of each particle 
is measured and the time to reach apogee is determined. The 
maximum height of a bounce is achieved when vertical veloc- 
ity is reduced to zero, which can be solved for the time to 2 o 
reach apogee. The total time of flight is twice that time, 
yielding the following solution for the initial velocity com- 
ponents in step 84 as a function of distance traveled and time 
of flight: 


Vyo = -a- 
Ax 

V X0 = — 


Equation (9) 


30 


Total velocity may be found by combining the vector com- 
ponents: 


v=Vv x 2 +v y 2 Equation (10) 

Ideally, the velocity decreases after each bounce by the 
COR of the particle and surface so that the reflected velocity 
in step 86 is: 

40 

=^COR partlcle COR surface Equation (1 1) 

Real particles, however, have angular momentum; energy 
transferred between translational and angular momentum 
will cause a reduction or increase in the velocity. Assuming 
the particle and surface COR are fixed properties, the ratio of 45 
velocities at each bounce (v^Vq, v 2 /v 15 etc) will be a constant 
plus or minus the fraction of energy participating in the trans- 
lational — angular momentum conversion in step 90. The 
energy balance is therefore: 


1 7 

K + U = -mv + Ioj 

2 

\ mv \ + Ioj\ = j mvl + I a)l - i m[(COR)v 0 ] 2 


Equation (12) 


55 


From which we can get the predicted velocity v x or mea- 
sured change in rotational energy: 


60 


Vl 



AIto 2 = -m{[\-{CORf]vl-v\) 


Equation (13) 


65 
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If the changes in rotational energy are distributed normally, 
the Coefficient of Momentum transfers (step 100) may be 
considered the standard deviation of the energy change nor- 
malized by kinetic energy: 


COM = stdevl 


A lof 
^~~ 2 


Equation (14) 


Furthermore, assuming the changes in angular velocity are 
distributed normally, the particle COR may be found by using 
equation (4) as the average of many measurements. The indi- 
vidual COR values calculated this way will be normally dis- 
tributed around the particle COR due to the effect of the 
momentum transfer. Then, using a surface for which 
COR SUK f ace is known (such as 0.9 for a glass plate), equation 4 
may be easily solved to determined COR particle as shown in 
step 102. 

Referring now to FIG. 6, the methodology for obtaining the 
diffuseness parameter will be described. The methodology 
begins in step 104 by bouncing a plurality of particles (one at 
a time) on a plate under the influence of gravity and the 
location and times of impact are observed, for example with 
a pair or high speed cameras. 

The angle of incidence when impacting after the ballistic 
trajectory it followed equals the launch angle. The variation in 
the launch angle after impact is modeled according to the 
following relationship using diffuseness to interpolate 
between the specular direction and a random lambertian vec- 
tor: 


V ! = v^^(l-diffuseness)+ v /aiM ^ fiaM (diffuseness) 


0=arc cos( v v Equation (15) 

The velocity and direction of the incoming particle is deter- 
mined (step 106) from the measured data and then the veloc- 
ity and direction of the reflected particle is determined (step 
108) from the measured data. Next, in step 110, the angle 
variation may be found from the measured data using: 


a ( v xiVxo + vyiVyo \ Equation (16) 

u = arccosl — ( t I 

V V V il + V Yl V V io + V Y0 ) 

Since the maximum deviation would occur when the lam- 
bertian vector is parallel but opposite the X velocity compo- 
nent, and the specular direction has the initial velocity com- 
ponents, 


vector = <vxo,vro) Equation (17) 

Vl = <v* 0 , Vyo)(l -d) + <-v*o, 0 )(d) 

= <(1 - 2 d)v X 0 , (1 - d)v Y o) 

= <V X l, v n ) 

d =i(i-tn) = (i_ni) 

2 V Vxo' ' Vyo 

Diffuseness can then be solved for using either horizontal 
or vertical velocity change, assuming enough trials have been 
run so that the worst case deviation from specular has been 
measured. 
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Alternatively, diffiiseness can be solved for by geometric 
reasoning. 

For small particles, it is possible to find sticking velocity by 
reducing the drop distance until no bounces occur and calcu- 
lating the velocity from the drop distance. Alternatively, the 
model may implement a velocity dependent sticking velocity 
as proposed in the Thornton et al. paper discussed above. 

One application of the present invention is spacecraft con- 
tamination control. Particle contamination can cause failure 
of spacecraft mechanisms and degradation of optical imaging 
capability, and can interfere with guidance and navigation 
sensors. By implementing the particle surface interaction 
model of the present invention, one can predict contamination 
transport, or where particle contamination will likely 
develop. By knowing the relative movement of the spacecraft, 
rotational and translation movement in orbit, the likely 
sources of particle contamination, the interaction model can 
be integrated into an algorithm to map particle contamination 
on the surfaces of the spacecraft. Knowledge of likely con- 
taminated surfaces is useful during the design and control of 
spacecraft. By continuing to apply the interaction model, the 
design can be altered and evolve to curb or isolate contami- 
nation on certain surfaces of the spacecraft. Particle surface 
contamination can also be a factor to control flight dynamics 
of an orbiting space, again curbing or isolating particle con- 
tamination on more sensitive portions of the spacecraft. 

It should be understood that while the bulk of the foregoing 
discussion has focused on particulate contamination on 
spacecraft, the present invention is not limited to such appli- 
cations and is applicable to other particle sources and envi- 
ronments. 

FIG. 7 is a block diagram that illustrates an exemplary 
hardware environment of the present invention. The present 
invention is typically implemented using a computer 114 
comprised of microprocessor means, random access memory 
(RAM), read-only memory (ROM) and other components. 
The computer may be a personal computer, mainframe com- 
puter or other computing device. Resident in the computer 
114, or peripheral to it, will be a storage device 116 of some 
type such as a hard disk drive, floppy disk drive, CD-ROM 
drive, tape drive or other storage device. 

Generally speaking, the software implementation of the 
present invention, program 118 in FIG. 7, is tangibly embod- 
ied in a computer-readable medium such as one of the storage 
devices 116 mentioned above. The program 118 comprises 
instructions which, when read and executed by the micropro- 
cessor of the computer 114 causes the computer 114 to per- 
form the steps necessary to execute the steps or elements of 
the present invention. 

It will be apparent to those skilled in the art having regard 
to this disclosure that other modifications of this invention 
beyond those embodiments specifically described here may 
be made without departing from the spirit of the invention. 
Accordingly, such modifications are considered within the 
scope of the invention as limited solely by the appended 
claims. 
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The invention claimed is: 

1. A program storage device readable by a machine, tangi- 
bly embodying a program of instructions executable by the 
machine to perform method steps for determining particle — 
5 surface interactions, the method steps comprising the steps 
of: 

inputting particle and surface properties; 

tracking a plurality of particles through capture on a sur- 
face; and 

determining a distribution of particles reflecting from the 
surface; 

wherein the step of tracking a plurality of particles 
through capture on a surface comprises the step of 
initializing each particle to a position and velocity and 
moving the particles and interacting with a surface; 
15 wherein moving the particles and interacting with a 

surface comprises the steps of 
determining a path of each of the particles; 
checking whether each of the particles intersects with 
the surface and proceeding to the next step for each 
20 particle where there is an intersection; 

checking whether each of the particles is captured by 
the surface and proceeding to the next step for each 
particle which is not captured by the surface; and 
determining a new direction and velocity for each 
25 particle; 

wherein the step of determining a new direction and 
velocity for each particle comprises the steps of: 
determining a coefficient of restitution; 
determining a momentum loss from the coefficient of 
30 restitution; 

determining a pre-impact kinetic energy; 
determining a post-impact kinetic energy; 
determining a coefficient of momentum transfer; 
determining a change in rotational momentum of each 
35 particle; 

updating the post-impact kinetic energy and rota- 
tional momentum of each particle; 
determining a post-impact velocity of each particle; 
determining a total momentum of each particle; 

40 determining a diffuseness of each particle; 

generating a diffuseness vector for each particle; 

determining a specular reflection vector; 

determining a post-impact reflection direction; and 
determining a reflection velocity vector 
45 wherein the step of determining a diffuseness of each 
particle comprises the steps of: 

inputting measurements of an apogee of each bounce of a 
plurality of particles bounded on a plate and time to 
reach the apogee; 

50 determining a velocity and a direction of each incoming 

particle; 

determining a velocity and a direction of each reflected 
particle; 

determining an angle between a specular reflection vec- 
55 tor and an actual reflection vector; and 

determining diffuseness from the angle. 



